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In these proceedings, I shall review the basic concepts of perturbative QCD in its high-energy 
limit. I shall concentrate on the approach to the unitarity limit, usually referred to as saturation, as 
well as on the gluon-number fluctuations the importance of which has recently been discovered. I 
shall explain the basic framework showing the need for those phenomena, first, from a simple picture 
of the high-energy behaviour, then, giving a short derivation of the equation driving this evolution. 
In the second part, I shall exhibit an analogy with statistical physics and show how this allows to 
derive geometric scaling in QCD with saturation. I shall finally consider the effects of gluon-number 
fluctuations on this picture and draw the physical consequences, i.e. a new scaling law, arising from 
those results. 
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I. INTRODUCTION 

The quest for the high-energy behaviour of perturbative QCD started thirty years ago, soon after QCD was proposed 
as the fundamental theory of strong interactions. To grasp this problem, it has been realised that saturation effects 
were to be considered in order to satisfy unitarity constraints. Throughout these proceedings, I shall summarise the 
present status of our understanding of this limit, which basically means finding an equation giving the evolution of 
amplitudes towards high energy and finding the properties of its solutions. 




log(Q 2 ) 

FIG. 1: Picture of the proton in DIS 

In order to motivate the physical picture we want to reproduce, let us consider the problem of Deep Inelastic 
Scattering (DIS) i.e. -y*p — > X collisions. This process depends on two kinematical variables: the virtuality Q 2 of the 
photon, imposing l/Q as the (transverse) resolution scale at which the photon scans the proton, and the Bjorken x, 
given by the fraction of the proton momentum carried by the parton struck by the virtual photon in a frame where 
the proton is moving fast. This last variable is related to the centre-of-mass energy s = Q 2 /x, meaning that the 
high-energy limit corresponds to the small- x limit (at fixed Q 2 ). For practical purposes, we also introduce the rapidity 
defined through Y — log(l/x). 

In Figure ^ we have represented the typical configuration of the proton in different domains of this phase space. 
We start with a proton at low Q 2 and low energy, represented as three partons (bottom- left part of the picture). If 
we increase the virtuality of the photon, we are able to resolve its partons into smaller ones. This type of evolution 
is described in perturbative QCD by a renormalisation group equation known as the Dokshitzcr-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equation In that limit, the number of partons rises logarithmically while their typical 
size decreases like l/Q 2 so that the proton becomes more and more dilute. In what follows, we shall be interested 
in the other type of evolution, namely the evolution towards smaller values of x for constant Q 2 . If one boosts the 
proton, what typically happens is that we create new partons of a size comparable to the parents' ones. Obviously, 
during this evolution, known as the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution 0, the proton approaches the 
black-disc limit for which partons start to overlap. At that level, they start to interact among themselves and unitarity 
corrections are to be considered 3J. This blackening of the proton is called saturation and we shall describe it in 
details from perturbative QCD in these proceedings. 

In addition, it has recently been discovered that gluon-number fluctuations play a crucial role in the evolution 
towards high energy. As we shall see later on, those fluctuations, dominating the evolution in the dilute tail (i.e. at 
large momentum) of the system, are amplified by BFKL evolution and thus influence greatly the whole evolution. 

In the first part of these proceedings, we shall derive the equations describing the evolution towards high-energy. 
We shall begin our analysis with a short introduction to the dipole picture which provides a very intuitive framework 
to deal with high-energy QCD at large-A^. This will allow for a derivation of the BFKL equation. The solution 
of the BFKL equation naturally leads to the introduction of unitarity corrections and to new evolution equations - 
the Balitsky/JIMWLK hierarchy 0,0 and the Balitsky-Kovchegov (BK) equation We shall then derive and 
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discuss the contribution due to gluon- number fluctuations [3, la, ul LL0| • 

In the second part, we shall consider the solutions to those evolution equations. First, we shall show that, within rele- 
vant approximations, the BK equation reduces to the Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation 
[l2|. well studied in statistical physics. This allows for a description jlH of the asymptotic behaviour of the scat- 
tering amplitudes in terms of travelling waves, translating in QCD into the geometric scaling property |l3j . We 
shall also show that geometric scaling extends to nonzero momentum transfer 0, . In the last section, we shall 
come back to the effect of gluon-number fluctuations, which have recently proven to have important consequences 
[Til IvA fljl liH I20L I2TL I22I l23j| on the approach to saturation. 



II. EVOLUTION TOWARDS HIGH ENERGY 



A. Small-x gluons and the dipole picture 



In order to clearly emphasise the need for a resummation at high energy, let us start with a fast-moving quark of 
momentum p. This quark emits Bremsstrahlung gluons characterised by their transverse momentum and their 
longitudinal momentum k z — xp. The probability for emitting such a gluon is 



dP oc a 



dx dk^ 

~"ky 



When we consider large transverse momentum (high Q 2 and fixed x in DIS), the collinear divergence dk^/k^ needs 
to be resummed and this gives rise to a renormalisation group equation known as the DGLAP equation. Throughout 
this proceedings, we shall instead consider the limit of high energy, in which, we are sensitive to the gluons of small 
longitudinal momentum x <C 1 at fixed Q 2 . This leaves a large phase-space to allow for successive soft gluon emissions 
satisfying x <C x\ <C ■ ■ ■ <C x n <C 1 which gives a contribution of order 



■ f 1 — ... f — = -a n log(l/x)™. 



Even when the coupling a s is small enough to ensure applicability of perturbation theory, for sufficiently small values 
of x (i.e. at sufficiently high energy), we have a s log(l/a;) ~ 1 and all these contributions have to be resummed. 
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FIG. 2: Onium wavefunction 

In practice, it is more convenient to start with a quark-antiquark pair of transverse coordinates x and y. Successive 
soft gluon emissions can be considered as independent and are depicted in figureEl In order to simplify the discussion, 
we shall consider the large- iV c limit (N c is the number of colours). In that limit, a gluon of transverse coordinate z 
can be considered as a quark-antiquark pair at point z. This means that, instead of considering a wavefunction made 
of quarks, antiquarks and gluons, it is sufficient to consider an onium which is an evolving system of colourless qq 
dipoles. We can thus think only in terms of dipoles. This is the dipole picture introduced by Mueller 24] . 



FIG. 3: Gluon emission in the dipole picture is equivalent to dipole splitting 



At high-energy, all possible gluon emissions from one dipole are thus tantamount to this dipole splitting into two 
child dipoles as depicted in figure [21 We see that the probability density for a dipole of transverse coordinates (x, y) 



FIG. 4: BFKL evolution through dipole splitting in the projectile 



to split into two child dipole (x,z) and (z,y) can be computed from gluon emissions from the quark and antiquark 
lines. The computation is performed in old-fashioned light-cone perturbation theory but we shall simply give the final 
result 1 : 



dP = — M xvz dYd 2 z 



(1) 



with a — a s N c /Tr and 



y - z 
(y - z) 1 



(x-y) 2 

(x-z) 2 (z-y) 



(2) 



B. The BFKL equation 

Let us denotes by (T xy ) the scattering amplitude 2 for a dipole made of a quark at transverse coordinate x and an 
antiquark at y. Getting the evolution of (T) from the dipole picture is pretty straightforward. Indeed, if one boosts 
the dipole (x, y) from a rapidity Y to a rapidity Y + SY, this dipole shall split into two dipoles (x, z) and (z, y). Each 
of those dipoles can interact with the target which leads to the following evolution equation 3 

dy (T xy ) = |- J M xyz «T„) + (T zy ) - (T xy )) (3) 

where we have explicitly used the probability density for dipole splitting. Equation |j3J is the BFKL equation, derived 
in the mid-seventies Q by Balitsky, Fadin, Kuraev and Lipatov. 

Note that here we have derived the equation by putting the evolution in the projectile dipole. It is also possible to 
write the evolution equation for the dipole density n xy in the target and to obtain |J3) by letting the projectile dipole 
interact with the evolved target. We sketch out the derivation of the BFKL equation for the dipole wavefunction of 
the target in Appendix [X] 

This equation being linear, one expects its solution to grow exponentially. To be more explicit, if one neglects 
impact-parameter dependence, i.e. assumes that (T xy ) = (T(r = |x — y|)), the solution of © is found by Mellin 
transform: 

(T(r))= /ilroW^W^ 10 ^/^, (4) 

J 2l7T 

where x{l) — 2^(1) — ^(7) — "0(1 — 7) is the BFKL eigenvalues plotted on Figure 03 and To (7) describes the initial 
condition. For a fixed dipole size, one finds the high-energy behaviour by expansion around the saddle point 7 = 1/2. 
One obtains 



(T(r)) ~ - exp 
ro 



2 x "(i/2)aY_ 

with cop = 47rlog(2)a being the BFKL Pomeron intercept. This expression suffers from two major problems: 



1 The first equality shows explicitly the two contributions coming from emission by the quark and the antiquark. 

2 The notation (•) denotes the average over all possible realisations of the target colour field. It can be understood as the expectation 
value for the T-matrix operator over the target wavefunction. 

3 The last term, coming with a minus sign, corresponds to virtual corrections 
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• first, even if one start with an initial condition peaked around a small value of r, when rapidity increases, the 
amplitude starts to diffuse to the large dipole sizes i.e. to the non-perturbative domain; 

• second, the solution of the BFKL equation grows exponentially with rapidity. Hence it violates the unitarity 
constraint T < 1 obtained from first principles. 




2 - 

I 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

7 

FIG. 5: The BFKL eigenvalues 



C. Saturation and the BK equation 

This unitarity problem is precisely the point where we meet the requirement to take into account multiple inter- 
actions. The most straightforward way to see this is to come back to the splitting of a dipole (x, y) into (x, z) and 
(z, y). When the target becomes dense enough, both dipoles can scatter on it. We thus have to take into account 
this contribution, represented in figure^] This leads to a quadratic suppression term in the evolution equation which 
becomes 

dr < T xy> = £ J M xyz (<T„) + (T zy ) - (T xy ) - (Ti 2 z ) zy )) . (5) 

The first thing to remark is that this new term is of the same order as the previous ones when T 2 ~ T or, equivalently, 
when T ~ 1. It is thus, as expected, a mandatory contribution near the unitarity limit. However, equation (JSj involves 
a new object, namely (T xz . zy ) which probes correlations inside of the target. In a general framework, one should then 
write down an equation for (7 12 ), which will involve {T 2 ) through BFKL-like contributions and (T 3 ) from unitarity 
requirements. This ends up with a complete hierarchy, giving the evolution for each (T k ^, known as the (large-iV c ) 
Balitsky hierarchy Q. 

If the target is sufficiently large and homogeneous, one can simply assume (T 2 z;zy ) = (T xz ) {T zy ), ending up with 
a closed equation for (T) 

dy ( T * y ) = £ / -Mxyz ((Txz) + (T zy ) - (T xy ) - (T„) (T zy » . (6) 





FIG. 6: Multiple scattering between a target and a projectile made of two dipoles. 
in the target. 



This can be interpreted as Pomeron merging 
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log(fc 2 /^) 

FIG. 7: Dipole occupation number shown as a function of the dipole momentum k for Y = 0,4, 8, 12, 16, 20. The dashed lines 
correspond to the solution of the BFKL equation while the solid lines correspond to BK evolution. 




FIG. 8: Diagrams contributing to the evolution of (T 2 ). They correspond, from left to right, to BFKL ladders, saturation 
corrections and fluctuations effects. 

This last expression is the Balitsky-Kovchegov (BK) equation 0, Q, which is the most simple equation one can 
obtain from perturbative QCD by including both the BFKL contributions at high-energy and the corrections from 
unitarity. Remembering that the dipole splitting is equivalent to gluon emission, we see from the diagram in figure 
El that the BK equation resums fan diagrams in addition to BFKL ladders. It is easy to check that (T) = is an 
unstable fix point of the BK equation, while (T) = 1 is a stable fix point, ensuring unitarity is satisfied. 

In addition, this new term also solves the problem of infrared diffusion. This can be seen on Figure [71 which shows 
the dipole occupation number for numerical solutions of the BFKL and BK equations in momentum space (k is the 
dipole momentum, canonically conjugated to its size). The BFKL solution exhibits the Gaussian shape predicted 
from (0J, extending both to the infrared and to the ultraviolet. In the BK solution, we clearly see that the infrared 
evolution is cut and the maximum of the distribution provides a natural scale Q S (Y), called the saturation momentum, 
increasing with rapidity. In general, one can say that saturation corrections cut the emission of dipoles of sizes larger 
that the inverse of the saturation momentum. We shall come back with more details to the study of the solution of 
the BK equation in the next section. 

D. Fluctuations 

Very recently, it has been realised that the question of the high-energy limit of QCD, even at the leading logarithmic 
accuracy, is not fully described by the Balitsky hierarchy. To see this more precisely, let us consider in more details the 
evolution equation for (T 2 ). It corresponds to the scattering of two dipoles off a target and the diagrams contributing 
to one step of rapidity evolution are shown in figure The first graph on the left of this figure is simply the BFKL- 
ladder contribution for which dy (T 2 ) oc (T 2 ^. The diagram in the middle corresponds to multiple scattering (see 
also figureEJ). It accounts for unitarity corrections and gives a contribution dy (T 2 ) oc (T 3 ). If those two graphs were 
the only relevant ones in the evolution, one would obtain the (large- N c ) Balitsky equation. Nonetheless, it has been 
argued that a third contribution, represented by the rightmost diagram of figure 03 has to be included. This graph 
describes multiple scattering off the projectile or, equivalently, gluon-number fluctuations inside the target. After the 
Pomeron-merging term, contributing to saturation, one thus now also includes the Pomeron-splitting contribution, 



FIG. 9: For a dilute target, the fluctuation diagram can be obtained from dipole splitting in the target. 



which gives rise to a new term in the evolution equation dy (T 2 ) oc (T). By combining Pomeron mergings and 
Pomeron splittings it is possible to build Pomeron loops. For that reason, in what follows we shall refer to those 
equations as the Pomeron-loop equations. 

For large- N c and at the two-gluon-exchange level, this new term has been computed We shall only give the 
main steps hereafter. We start with the observation that gluon-number fluctuations are expected to be important 
when the target is still dilute, which is always the case for sufficiently small dipole sizes. In that case, it can be 
understood as a superposition of dipoles and the s-channel gluon in the fluctuation diagram (last diagram of figure 
IHJ can be considered as emitted through a dipole splitting in the target, as represented on figure Since we work in 
the two-gluon-exchange approximation, the scattering amplitude and the dipole density are related through 



<T xy ) = a 2 s f .4 (xy|uv) 

J uv 



with Aq being the dipole-dipole amplitude at the two-gluon-exchange level 



A)(xy|uv) = - log 2 
This relation can be inverted to obtain 



(x-u) 2 (y- v) 2 



(x- v) 2 (y - u) 



(7) 



(8) 



-2 V 2 V 2 irp 
v x v v \^x 



The new contribution to the evolution of (T 2 iyi . x . 3y9 ) is thus 



8 / T 2 \i = ±JL(2!l\ 

t ^\- t x 1 yi;x aya /| fluot 2 27rV27ry 

/ M U vz^o(xiyi|uz)A(x2y 2 |zv)V 2 1 V 2 (T uv ) + (1 «-> 2), 

J UVZ 



(9) 



The correction © becomes of the same order as the BFKL contribution when T 2 ~ ct 2 T, or T ~ a 2 , i.e. in the 
dilute regime where, as expected, fluctuations should lead to important effects. Although this may at first sight seems 
irrelevant for the physics of saturation, one has to realise that, due to colour transparency, the amplitude becomes 
arbitrarily small for small dipole size, hence, even for dense targets, there is always some region of the phase-space 
where the system is dilute and its evolution governed by fluctuation effects. 

Let us give an additional argument in favour of the importance of the fluctuation term on high-energy evolution. 
Once we have a Pomeron splitting coming from a fluctuation, it grows, from BFKL evolution, like two Pomerons, 
i.e. T ~ a 2 exp(2utpY). This has to be compared with the one-Pomeron exchange T ~ exp(wpF). We see that 
at very large energies Y > ujp 1 log(l/a 2 ), BFKL evolution compensates the extra factor of a 2 coming from the 
initial fluctuation. We shall see more precisely in the last section that those fluctuation effects in the dilute tail have 
important consequence on the saturation physics. 

So far in this section, we have only discussed the evolution equation for {T 2 ). What we really get is an infinite 
hierarchy, where the evolution of (T) is given by (JSJ), as in the Balitsky hierarchy, and the evolution of (T k ) contains 
three contributions: a linear BFKL term proportional to (T k ) 1 a (negative) saturation term proportional to (T k+1 \ 
as it appears in the Balitsky hierarchy, and the new fluctuation contribution proportional to (T^ -1 ), similar to J^J. 
Remarkably enough, it is possible to show that this infinite hierarchy for the average amplitudes can be rewritten as 
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a single Langevin equation 4 

d Y T xy = £ J^M xyz [T xz +T zy -T xy -T xz T zy ] (10) 

where is a Gaussian white noise satisfying 

("uv Z; r) = 0, 

<^uvz;F^uW;y) = <KaF - aY')5™ ( U _ v ') 5 (2) (z _ z r )d (2) (v _ u / } _ 

Different realisations of the noise term lead to different events T Xi y, from which one can compute (T fc ) using standard 
statistical techniques. However, dealing with equation (|1(J|) remains very difficult since the noise term is off-diagonal 
and non-local. 

Finally, let us conclude this section by quoting that the equation (j^J), including linear, saturation and fluctuation 
effects is the most complete equation known in perturbative high-energy QCD so far. Its extension beyond the 
large- N c approximation is still a challenging problem. 



E. QCD as a reaction-diffusion process ? 

Let is make a short digression and consider a reaction- diffusion system in which we have a set of particles evolving 
through splitting at a rate 7 and merging at a rate a: 

A^A + A. 

a 

For simplicity, we shall not consider any spatial dimension, though it only introduces technical difficulties without 
altering the conclusions. If we want to study the time evolution of such a system, the best way is to start with 
the master equation which gives the evolution of P n , the probability to have n particles. It contains four terms, 
corresponding to gain and loss terms for splitting and merging weighted by adequate combinatorial factors: 



d t Pn = 7 (JV - 1)Pjv_i - 7 NP N + a N(N + 1)P N+1 - a N(N - 1)P, 



JY 



gain loss gain loss 

In practice, this equation is only helpful to find evolution of observable, averaged, quantities. The most simple 
example is the average particle number 



N=0 

which can be generalised to a subset of k particles 

Using the master equation, one finds the evolution of the particle number and its correlators after a bit of algebra 
d t (n k ) = 7 k (n k ) +^k(k- 1) {n^ 1 ) - a k{k + 1) <n fe+1 ) -ak (n k ) . 

More interesting is the evolution of the scattering amplitude for this system off a generic target, defined by 

00 

Ait) = (n k ) to (T k ) t _ tQ , 

fc=0 



4 We show the equivalence between the two formulations in appendix [E] for a simpler of Langevin equation encountered in Section 4. 
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where we evolve the particle system up to a time to and the scattering matrix for the target from to to t. To be 
consistent, this definition should not depend on any specific choice for t . This gives a renormalisation group equation 
from which one can infer the evolution of the scattering amplitude: 

dt (T k ) = 7 (T k ) ~ 7 (T fc+1 ) + a [(T^ 1 ) - (T k )] . 

BFKL sat. fluct. 

We clearly recognise in this expression the BFKL, saturation and fluctuations contributions obtained in QCD 5 . 

Inserting the correct spatial degrees of freedom to account for transverse coordinates, we find [25| that the high- 
energy QCD evolution hierarchy obtained in the previous section can be seen as an effective dipole reaction- diffusion 
process where rapidity plays the role of time and with the vertices 

7(xy^xz, zy) ~ aM xyz , (11) 
cr(xiyi,x 2 y 2 -► uv) - aa 2 s \7 2 u \/l [M uvz A {xiyi\uz)A {x 2 y 2 \zv)} (12) 

for splitting and merging respectively. One recognise in these expressions that the splitting is given by the BFKL 
kernel while the merging comes from the computation in the previous section. 

This interesting analogy is however to be taken with care. Indeed, the merging vertex is not positive defined. 
Physically, this comes from the fact that fluctuations really involve the gluonic degrees of freedom. Their interpretation 
as a reaction-diffusion system in terms of dipoles can hence only be effective. 

III. PROPERTIES OF THE BK SCATTERING AMPLITUDES 

In this section, we shall sketch out the main properties of the solutions of the evolution equations towards high 
energy. As we shall explain in details in the next lines the major source of information comes from the analogy 
between the BK equation and the Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation which is well studied 
in statistical physics. This section shall be devoted to derive this analogy and its consequences in the case of the BK 
equation. We shall first consider the impact-parameter-independent BK equation and then show how the arguments 
extend to the full BK equation. We leave the discussion concerning the effect of fluctuations for the next section. 

A. Statistical physics and geometric scaling 

So, let us start with the BK equation ©. We shall first restrict ourselves to the impact-parameter-independent 
version of the equation, for which one can easily go to momentum space by using 6 

1 f d 2 r P°° rlr 

T ( fc ) = V / T e T ( p ) = / - Jb(*r)T(r), 
2n J r z Jo r 

and where the equation becomes 

drf(k) = «( d 4 

7T J P Z 

It is actually possible to rewrite the integral in is equation in term of a differential operator: 

d Y f(k) = a X {-d L )f(k) - af 2 (k), (13) 

with L = log(fc 2 /fcg), fco being a soft reference scale and x(t) ^ s the BFKL kernel. 

In order to simplify the problem, we shall work in the saddle point approximation (often referred to as the diffusive 
approximation) i.e. expand the BFKL kernel to second order around 7 = 1/2, transforming the complicated differential 
operator in equation 1)13(1 into a second-order operator 

x(-dL)T(fc)«x(i) + lx"(!)09L + i) 2 . 



p 2 f (p) - k 2 f(k) k 2 f(k) 



\k 2 — p 2 



y/Ap* + fc 4 



af 2 (k). 



5 The second fluctuation term disappears in the case of QCD. 

6 In this section, we omit the (■} quotations as they do not play any role in the mean-field approximation. 



10 




FIG. 10: Formation of a travelling wave from a schematic evolution of a wavefront in time/rapidity. 



Up to a linear change of variable switching from Y and L to time t = aY and space x — L+cst.K and a renormalisation 
of the amplitude T — > u — cst.T, the BK equation in the saddle point approximation becomes 

d t u(x, t) = d^u(x, t) + u(x, t) — u 2 (x, t). 

This equation is nothing but the F-KPP equation studied in statistical physics since forty years and applying to 
many problem. It describes reaction diffusion processes in the mean-field approximation where one can have creation 
and annihilation of particles locally and diffusion to neighbouring site. This equation has many applications e.g. in 
chemistry and biology. 

One knows that the F-KPP equation admits travelling waves as asymptotic solutions i.e. there exists a critical 
velocity v c and a critical slope j c , determined only from the knowledge of the linear kernel of the F-KPP equation 
9^ + 1, such that (see figure ITT1 for a pictorial representation coming from the BK equation which we shall explain a 
bit later) 

u(x,t) ^ u(x - v c t) x> ~ at e-^ x - v ^(x - v e t). 

Although working in the diffusive approximation makes the discussion easier due to the direct matching with the 
F-KPP equation, this assumption is not required. Indeed, for an equation with a more complicated kernel (than 
d\ + 1), one can prove the existence of travelling waves provided the three following conditions are satisfied: 

1. the amplitude has as unstable fix point and 1 as stable fix point, 

2. the initial condition decreases faster than exp(— j c L) at large L (see below for the general definition of j c ), 

3. the equation obtained by neglecting the nonlinear terms admits superposition of waves as solution: 

f ^7 

u { x M\in = i 77- Ml) exp [—y(x - v y t)] , 
J LlK 

where the wave of slope 7 travels at a speed w 7 . 

Then, travelling waves are formed during evolution with the critical parameters 7 C and v c obtained, from the linear 
kernel only, through the relation 

« c =x'(7c) = ^- 
7c 

This corresponds to the selection of the wave with the minimal speed x(l)/l an d the point where group and phase 
velocities are equal. 

Before going any further, let us give a qualitative argument explaining how the minimal speed is selected. To 
simplify the discussion, we shall consider the initial front at the left of figure EI] At an initial time, it captures 
the main features of the spatial dependence: saturation (first condition here-above) followed by a decrease with the 
critical exponent j c and finally a steeper decrease (second condition here-above). After one step of time evolution, 
each of the decreasing fronts shall evolve with its own speed (those fronts are in a dilute regime so we can use the 
third condition). The result of this evolution is shown on Figure ITU1 and we clearly see that the region over which the 
critical slope, i.e. the one with the minimal speed, is growing during the evolution. This selection mechanism is very 
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FIG. 11: Numerical simulation of the rapidity-evolution of the BK equation. We start at Y = from the leftmost amplitude 
and evolve to higher Y using (|13|) . Amplitude is shown for Y = 5, 10, 15, 20, 25 and clearly exhibits a travelling-wave pattern. 




FIG. 12: Rapidity evolution of the saturation scale extracted from figure fTTI More precisely, the speed of the wave, i.e. the 
exponent of the saturation scale, is plotted. It goes to a constant as it should. 



general. It mainly relies on the properties of the linear equation and saturation only ensures that the slope is j c at 
some point in the initial condition. 

In the case of the BK equation, the first condition is satisfied due to BFKL growth and saturation, the second 
comes from colour transparency (T oc 1/fc 2 = exp(— L)) and the last one is equivalent to the solution Q (with r/ro 
replaced by ko / k) of the BFKL equation obtained by neglecting nonlinear terms in (|13[) . 

The critical parameters (for the complete BFKL kernel i.e. for (|13fl ). arc found to be 7 C 0.6275 and v c « 4.8836a. 
Written in terms of the QCD variables Y and k 2 , one can then write the asymptotic solution for the impact-parameter- 
independent BK equation under the form 



T{k;Y) Y =°° T 



/ e \ 


fe»Q s 


" k 2 - 


-7c 

log 


" k 2 ' 


\QI(Y)J 




lQ 2 s(Y)_ 


[Q 2 s(Y)\ 



(14) 



with the saturation scale given by 



Q 2 S (Y) Y ^ fc 2 exp 



v c Y-— log(Y) 
27c 



(15) 



The formation of a travelling-wave pattern when energy increases can easily be seen on numerical simulations of 
the BK equation. As displayed in figure ITT1 if one start with a steep enough initial condition (leftmost curve), the 
amplitude increases with energy and a wave moving into the dilute domain gets formed. From that simulation, one 
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can extract the saturation scale by solving T(Q 2 , Y) — N at each values of Y and for a fixed threshold N. In figure 
1121 we have plotted \og(Q 2 )/Y which goes to a constant value as it is expected from (|15fl . 

Equation (|14|l has a remarkable property: it proves that at high energy the amplitude, a priori a function of both 
Y and k, depends only on the ratio between k and the saturation momentum. This property, known as geometric 
seating, has been observed in the HERA measurements of the proton structure function. The fact that geometric 
scaling can be derived from the BK equation is one of the most important indication for the experimental observation 
of saturation. The saturation scale obtained from the structure function data is of the order of 1 GeV for x ~ 10 -5 . 
This means that high-energy QCD can indeed be studied from the point of view of perturbation theory. 



B. Nonzero momentum transfer 



Up to now, we have only discussed the impact-parameter-independent BK equation. We might therefore ask whether 
or not these arguments extend to the full equation including all phase-space degrees of freedom. A dipole of transverse 
coordinates (x, y) is then represented through its size r = x y and impact parameter b = (x + y)/2. The problem 
is then that dipole splitting is non-local in impact parameter. This means that the BK equation couples different 
values of b and we can not apply directly the previous arguments for each value of b. Again, the solution consists in 
moving to momentum space and replace the impact parameter by the momentum transfer q 



f(k,q)= / d 2 xd 2 ye *.x e ,(q-k). y 

J (x - y) 2 



The BK equation then takes a form [TJj for which the BFKL kernel is local in q (the non-locality of the non-linear 
term is not important for our purposes as its only role is to ensure saturation) 



drT(k,q) 



d 2 k' 



(k-k') 2 



?\k',q)- 3 



k /2 



f. / d 2 fc'r(k,k')r(k-k', q -k'). 



(q-k) 2 
(q-k') 2 



(16) 



T(k,q) 



We can now proceed with this equation in a similar way as for the the impact-parameter-independent case. The 
existence of travelling waves requires the three conditions stated previously to be satisfied. The two first ones are 
trivially satisfied 7 . For the third condition to be valid, we need to find solutions of the linear part of (jl6|) . i.e. the 
complete BFKL equation, which can be expressed as a superposition of waves. This consists in a careful treatment of 
the solutions of the BFKL equation 26] including all phase-space variables. It turns out that, using the momentum 
transfer q instead of the impact parameter b , leads to a powerful factorisation between the target and the projectile 
(see Appendix[C]for more details). Then, we can show |15J that when the dipole momentum k is much larger than the 
momentum transfer q and the typical scale of the target ko, the solutions of the BFKL equation are a superposition of 
waves and therefore, we obtain travelling waves for the full BK equation. More precisely, the high-energy behaviour 
of the amplitude takes the form 



f(k,q;Y) Y =°° T 



fc»Qa 



1 ~To 



Q 2 s (q;Y) 



log 



QUq;Y) 



This expression is exactly the same as for the previous case except that the saturation scale now depends on momentum 
transfer 



,(q;Y) = A exp 



v c Y 



27< 



log(F) 



7 Indeed, we have saturation and BFKL growth ensuring the first condition and colour transparency still apply for the second one. 
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with 




&0 if k > q 
q 2 if q > k ' 



We explain in Appendix El how the exponential behaviours required for travelling- waves formation appear and how 
the reference scale k or q emerges. Naively, we can expect the saturation scale at large q to be determined by q as, 
in the tail of the amplitude, we are dominated by the hard part of the scattering. In addition, when <;->0we should 
recover the ^-independent result and some soft reference scale should come back. 

It is very interesting to notice that the critical slope j c and speed v c , obtained from the BFKL kernel only, are the 
same as for the impact-parameter-independent case. 

The important point at this stage is that this study predicts geometric scaling at non-zero momentum transfer. 
Experimental measurements of the DVCS cross-sections or diffractive p-meson production are good candidates and 
can provide more evidence for saturation. 



IV. HIGH-ENERGY QCD WITH FLUCTUATIONS 

Let us now consider the effect of the fluctuation contribution © . This new term, added to the full Balitsky hierarchy, 
turns a simple non- linear equation into an infinite hierarchy with a complicated transverse-plane dependence 8 . We can 
however simplify the equation to a tractable problem by performing a local-fluctuations approximation. This amounts 
to simplify the dipole-dipole scattering amplitude Ao used to relate the dipole density n with the scattering amplitude 
T. We assume that two dipoles interact if they are of the same size and if their centre-of-mass are sufficiently close 
to allow for an overlap of the two dipoles. Within this approximation equation Q is replaced by 

(T xy ) = na 2 s |x-y| 4 (n xy ) , 

where k is an unknown fudge factor of order 1. 

Once this is done, the remaining steps are as follows: one Fourier-transform the dipole size r = |x — y| into 
momentum k and perform a coarse-graining approximation to get rid of the impact-parameter dependence of the 
fluctuation term. This computation results into the following simplified form of the hierarchy (we give only the 
expression for (T) and (T 2 ) for simplicity) 

d Y (T k ) = a X (-d L )(T k )-a(Tl k ), 
+ a na 2 k 2 S(k 2 — k 2 ) {T kl ) , 

where, as previously, Li — log(k? 



Again, this infinite hierarchy has the advantage that it can be rewritten under the form of a Langevin equation' 



d Y T(L) = a [ X (-d L )T(L) - T 2 (L) + ^Jkc^T{L) v {L, Y)\ , (17) 

where rj is a Gaussian white noise satisfying the following commutation relations: 

(V(L,Y)) = (18) 

(TjtLi.yiMia,^)) = -S(L 1 -L 2 )S(Y 1 -Y 2 ). 
a 

The Langevin equation is a stochastic equation describing an event-by-event picture. Each realisation of the noise 
term corresponds to a particular evolution of the target and one can show (see Appendix that once we average 
over those realisations using the correlations (|18f) the complete hierarchy for the evolution of (T k } is recovered. It is 
interesting to notice that equation (|17fl is formally equivalent to the BK equation with an additional noise term. 



The fluctuation term involves an awkward integration. By integration by part one can express it as a vertex applied to T but the 
integration over the internal variable z in J5J is unknown. 

As we have seen previously, the complete hierarchy using J^l can a l so be rewritten as a Langevin equation 1101 . However, the noise 
term here, being local, appears to be much simpler than for the general case 1101 . 
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FIG. 13: Speed of the wave for different values of the noise strength /taj. The BK speed has been added for comparison. One 
sees that the speed decreases when KO? a increases. 




FIG. 14: Numerical simulations of equation 1171 . The dashed curves are different events corresponding to the same initial 
condition with different realisations of the noise term. Travelling waves are observed for each curve, together with dispersion. 
The black curve is the average amplitude. 



As for the case of the BK equation, we can first restrict ourselves to the diffusive approximation. This leads to 
the stochastic FKPP (sFKPP) equation which amounts to add a noise term to the FKPP equation. It has found 
many applications in various fields, e.g. in the description of reaction-diffusion systems. The noise term appears once 
we have to consider discreteness effects {e.g. for a reaction-diffusion system on a lattice, only an integer number of 
particles are admitted per site). When the number of particles involved goes to infinity the mean-field approximation 
seems justified and the (non-stochastic) FKPP equation is valid. 

The basic effects of the additional noise term in the sFKPP equation are known, at least qualitatively. Despite 
the fact that the fluctuations are only expected to have large consequences on the dilute tail of the wavefront, these 
modifications change the picture at saturation. In order to test the validity of those results in the case beyond the 
diffusive approximation, i.e. with the full BFKL kernel, we have performed 0] numerical studies (see also ^tJ) °f 
the QCD equation 1)170. Moreover, many analytical results are only known in the limit where the noise strength na 2 s 
is (irrealistically for QCD) small while we have concentrated our numerical work of physically acceptable values. In 
the following paragraphs, we review the main effects of the fluctuation contribution and show that they are observed 
in the numerical analysis. 

First of all, for a single event, the evolved amplitude shows a travelling- wave pattern (up to small fluctuations 
in the far tail which, at least a this level, are irrelevant for the wave pattern at saturation). This means that each 
single realisation of the noise leads to geometric scaling, as it is the case for the BK equation. At this level, the 
major difference comes from a decrease of the speed of the wave. Analytically, this can only be computed in the limit 
KQj i.e. when the fluctuations are a small perturbation around the mean-field behaviour 10 . It has been found 



A recent analysis has also computed its strong- noise behaviour llSl . 
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FIG. 15: Dispersion (squared) of the position of the events as a function of rapidity. As expected, it dispersion increases like 
vF but few dispersion is obtained in early stages of the evolution. 
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FIG. 16: Evolution of the averaged amplitude for different values of the noise strength, compared with the mean field result 
(dotted curve). The curves correspond, from left to right, to Y = 0, 10, 20, 30, 40 and 50. 



[Hill that 

a7T 2 7 c x"(7c) 

v — > v c 5 . 

oJk^O 2 log (a 2 K) 

Unfortunately, this expression gives only reliable results for extremely small values of rea 2 — > such as 10 -20 , which 
is not sufficient for realistic situations. The numerical simulations we have performed clearly exhibits this decrease of 
the speed when the noise strength na 2 increases (see figure EU- 

The second noticeable effect of the noise term is to introduce dispersion between different events. Different events 
have the same shape but their position X(Y) (\og[Q 2 {Y)] in physical variables) fluctuates. At a given rapidity, one 
can compute the dispersion of those events. This diffusion process being very similar to a random walk, one expects 
that the dispersion of the events behaves like yY : 

(X 2 ) y - (X) 2 Y y ~°° D diS Y. 

The parametric dependence of the diffusion coefficient Ddiff has been computed numerically ljj and behaves 11 like 
| log _3 («;a 2 )| for small na 2 s . 

Concerning the numerical simulations, the dispersion of the wavefront between different events and its increase with 
rapidity are manifest on figure ITU One can then extract the dispersion as a function of rapidity. The resulting curve, 
shown on figure 1151 for different values of the noise strength, calls for two remarks. First, the expected asymptotic 



11 Recently, this coefficient has been computed analytically l'2(l for rect^ <C 1. 
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behaviour (dispersion oc \Y) is observed and the diffusion coefficient -Ddiff increases with rapidity and with the noise 
strength. However, for small values of the rapidity, we do not observe a significant dispersion. This lack of dispersion 
may come from the fact that we have to wait for the travelling front to get formed before dispersion becomes significant. 
However, those effects in the beginning of the evolution are perhaps not universal and deserve more detailed studies. 

This dispersion of the events has an important physical consequence: although each single event displays geometric 
scaling, once we compute the average amplitude, dispersion induces geometric scaling violations. This effect is best 
seen on figure IT^l where we have compared the evolution including the fluctuations to the BK results. While in the 
BK equation, a fixed travelling-wave pattern is formed, once fluctuations are taken into account, we see a broadening 
of the average amplitude as rapidity increases. 

Those violations of geometric scaling can be studied |2lj | in a model which encompasses all the physical arguments 
presented here-above and which shall prove to be very efficient for phenomenological studies. In order to reconstruct 
physical amplitudes, we need two building ingredients: 

• a front-by-front amplitude which displays geometric scaling. This amplitude has to satisfy unitarity in the 
infrared and must decrease exponentially in the ultraviolet. The simplest choice, which is actually sufficient for 
what follows, is 



Invent (P _ Pa) 



exp[-7o0- ps)] Hp>Ps, 
1 if p < p a , 



(19) 



where we have introduced the notations p = log(rp/r 2 ) and p s = log[l/(r 2 Q 2 )]. This expression shows indeed 
geometric scaling since it is a function of p — p s only. The precise value of 70 is not relevant for what follows, 
although it might be adequate to adopt 70 = 7c — 0.6275 or 70 = 1. 

• a dispersion between different events. Physically, this dispersion of the position of the events, i.e. of their 
saturation scale, means that the saturation scale becomes a random variable. We thus need a probability 
distribution for p s . We shall take a Gaussian form 12 



PiP.) 



V2 



exp 



(Ps - Psf 



2a 2 



(20) 



The parameters of the Gaussian (the average saturation scale p s and the dispersion a 2 ) are expected to grow 
linearly with rapidity, but we shall treat them as parameters. 



The average amplitude can be constructed from (|19J) and (|2()|> 

/oo 
dp s P{ps) T cvcnt (p - p s ) 
-oo 



= -erfc 
2 



V2a 



1 erfc 

2 



700- 



2 \ i 



exp 



^,2^.2 



with z = p — p s - In the last line, the first term comes from the event-by-event fronts which are at saturation (T = 1) 
while the second one comes from the exponential tail. Practically, there are two interesting limits to analyse in more 
details: 

1. 7pg 2 < 1 

When the dispersion is small (compared to 1 /70 which is the typical decay length of the exponential tail) , the 
dispersion can be neglected. We thus recover 



( T (p)) = T event (p- p s ), 



(21) 



and geometric scaling is satisfied. Since a 2 oc Y, this is what happens in early stages of the evolution. 



12 We can actually prove 122 
domain. 



that, at high-energy, the probability distribution can be considered as Gaussian over an arbitrarily large 
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2. 7qQ- 2 > 1 

For large dispersion, i.e. at very high energies, we find 

™> = r rfc (7fe)- < 22 > 

This expression is valid for —70c 2 < z < To^ 2 which is arbitrarily large at high energy. 

From this analysis, we deduce two very important physical consequences for the high-energy behaviour of the 
amplitudes: 

• a new scaling law, called diffusive scaling, appears. The amplitude scales like z/cr, which, in more physical 
variables means that it only depends on the ratio log(r 2 Q 2 ) / 'VY '. It is possible to show that this property of 
diffusive scaling extends to physical observables such that DIS, diffractive DIS 0,113 and gluon production 

mm 

• The error function in (I22|) fully comes from the event-by-event fronts which are at saturation. In other words, 
the high-energy behaviour of the amplitudes can be obtained by taking T C vent = 9{Ps — p)- Thus, at high energy, 
the scattering amplitudes are dominated by black spots, T = or T = 1. 

V. DISCUSSION AND PERSPECTIVES 

Let us now summarise the main points raised in these proceedings. As this document is in itself some kind of 
a summary, I shall only pick up the most important points, refer to various references for detailed approaches and 
discuss open questions. 

First, we have shown that it is possible to describe the evolution to high energy in pQCD by an infinite hierarchy 
of equations (see (01 )■ This gives the evolution of the (T) matrix (averaged over the target wave- function) and its 
higher-order correlations (T k y In this hierarchy, the evolution of ^T fe ) contains three types of contribution: 

1. the linear BFKL growth, proportional to a(T fc ). This is the usual high-energy contribution computed thirty 
years ago and corresponding to the exchange of k pQCD Pomerons. It leads to a fast (exponential) increase of 
the scattering amplitude. 

2. the saturation corrections, proportional to a (T k+1 }. This negative term becomes important when the amplitude 
reaches unitarity and it allows for the constraint T(r, b) < 1 to be satisfied. 

3. the fluctuations term, proportional to aa 2 (T k ~ 1 'j. These fluctuations correspond to gluon-number fluctuations 
in the target. They play an important role when the amplitude is of order a 2 or, equivalently, when the dipole 
density is of order one, i.e. in the dilute regime where one indeed expects fluctuation effects to appear. 

If one neglects the fluctuation contributions, we recover the Balitsky hierarchy (in its large- N c formulation) and, 
in the mean field approximation, the BK equation is obtained. This last equation, much simpler than the infinite 
hierarchy captures many important points concerning the physics of saturation. We also have to point out that, as 
long as we do not take into account fluctuations, the evolution at all orders in 1/N C exists under the form of the 
Balitsky/JIMWLK equation. The Balitsky equations take the form of a hierarchy, similar to the one introduced 
in these proceedings but involving, in addition to dipoles, more complicated objects like quadrupoles, sextupoles, 
... The JIMWLK equation is an equivalent formulation based on the Colour Glass Condensate approach, in which 
the probability density for the colour field of the target evolves through a functional equation. Further information 
concerning the Balitsky (resp. CGC) approach can be found in [2^ (resp. HjJ). 

Many additional things can be said concerning those different contributions, e.g. concerning the equivalence between 
descriptions from the target and projectile point of view (see e.g. H3), or the attempts to describe the evolution 
beyond its large— N c limit |2l|. Those questions are active fields of research and the interested reader is forwarded to 
the corresponding references for further information. 

In the second part of this overview of high-energy QCD, we have concentrated ourselves on the physical consequences 
arising from saturation and fluctuation effects. We have investigated how those contributions manifest themselves on 
the scattering amplitude and change the exponential behaviour obtained from the BFKL evolution. Again, the most 
important points can be summarised in two important steps which have their respective physical consequences: 
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1. Considering BFKL with saturation effects (mean- field picture), one can show that the BK equation lies in the 
same universality class as the FKPP equation. This implies that travelling waves are formed during the evolution 
towards high energy. From the BFKL kernel one obtains the critical parameters which are the anomalous 
dimension observed in the large-Q 2 tail and the speed of the wave, equivalent to the exponent of the saturation 
scale. In physical terms, these travelling waves correspond to geometric scaling. The experimental observation 
of geometric scaling in DIS at HERA can thus be seen as a consequence of saturation. It is extremely important 
to realise that geometric scaling is a prediction from saturation physics which extends to large values of Q 2 , 
far beyond the saturation scale itself 13 . This geometric scaling has also been predicted at nonzero momentum 
transfer, a result which may be applied for example to vector-meson production or DVCS. 

2. If one includes the effects of fluctuations, the evolution becomes a Langevin equation equivalent to the stochastic 
FKPP equation, including a noise term. For each realisation of this noise, we observe geometric scaling with a 
decrease of the speed w.r.t. the mean-field case. If we consider a bunch of events, we observe that, although 
they all have the same shape, their position is diffused with diffusion increasing like y/Y. This dispersion of the 
saturation scale for a set of events implies that geometric scaling is violated for the averaged amplitude. Recent 
studies indicate that geometric scaling is still valid at small rapidities. At higher energies, a new (diffusive) 
scaling is predicted [2lJ in which the amplitude scales w.r.t. the scaling variable r = \og[k 2 / Q 2 (Y)]/ \/Y . 

Although these two pictures endow most of the physical effects, they are far from being completely understood. 
Among the improvements, one can quote phenomenological applications like geometric and diffusive scaling predictions 
for vector-meson productions, for diffraction |2l| and for LHC physics |23j. On more theoretical grounds most of the 
work that still need to be done concerns the effects of fluctuations. Most of the results known so far are derived 
within the local approximation for the noise term and neglecting the impact parameter dependence. In addition, 
apart from numerical studies, analytical results can only be applied to irrealistically small values of a s and a better 
understanding for more physical values is still lacking. Although the existing picture is believed to contain all the 
qualitative effects, a more precise treatment is certainly an interesting and challenging problem to address. 

Within the framework presented here, the link between the QCD evolution equations and equations from statistical 
physics has led to a large number of interesting ideas and results. Even in statistical physics, many questions are yet 
opened, hence, extending the QCD picture including more sophisticated treatments of the BFKL kernel and more 
precise forms of the noise term is certainly not a straightforward task. The questions of making clear predictions for 
phenomenology, of understanding the effects of fluctuations beyond the local-noise approximation and of including 
impact-parameter dependence are hence expected to give interesting work in the near future. Among particle physics, 
high-energy QCD is thus one of the most active fields. Because of those open questions and because of the requirement 
for an excellent knowledge of QCD at the LHC, we expect that high-energy QCD is going to be even more active in 
the near future. 
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Let us show that it is also possible to obtain the BFKL equation by looking at the evolution of the dipole densities 
in the target. We hence need to compute how, within a small increase of rapidity, the density of dipoles of coordinates 
(x, y) evolves. The answer is 



The three contribution in this equation has well determined origin: the first term correspond to the situation where, 
at rapidity Y, we had a dipole (x, z) which splits into (x, y) and (y, z) hence the positive contribution to the evolution 
of n xy . The second term has similar explanation with an original dipole (z,y) and the third one corresponds to 



The geometric scaling window grows, in logarithmic units, like y/Y beyond the saturation scale. 
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the situation where we had one dipole (x, y) at rapidity Y which disappears through splitting into (x, z) and (z, y), 
leading to a loss term coming with a minus sign. 

Getting to the the evolution of the scattering amplitude when we probe that system with a dipole (u, v) required a 
few technical manipulations. First, we shall use eq. Q to construct the scattering amplitude from the dipole density 
by convolution with the fundamental dipole-dipole interaction 14 : 



dy <T UV ) = f -4o(uv|xy) [M 

J xyz 



xzy^x,z ~\~ >^zyx^z,y -^xzy^x,z] 



The trick to recover © is to use conformal properties. For example, for the first term in the equation, we perform 
the following change of variables 15 for the integration over y: 

, ay + b \a = —d = xu~zv, c = u — v+x — z, 

y^y= f{y) = — —7 with < 

cy + a I b = uvz — uvx + xzv — xzu, 

which has the effect to interchange x with u and z with v. It thus implies 

A(uv|xy).M xzy -> / ^o(/(u)/(v)|/(x)/(y))M /(x)/(z)/(y) 



y 



■Ao(xz|uy')M U vy 

y' 

Carrying on with the x and z integration will give 

/ M uvy ,(T uy ,). 
Jy' 

Renaming y' into z and performing a similar treatment for the remaining terms, one recovers the BFKL equation |J3j. 
APPENDIX B: LANGEVIN EQUATION AND INFINITE HIERARCHY 



In this appendix, we show that a Langevin equation with a noise term proportional to ^Ju(x, t)rj{x, t) (with 77 a 
Gaussian white noise), turns into the fluctuation term in the hierarchy for averaged amplitudes. To stay as close as 
possible to the problem we are concerned with, let us consider a sFKPP-like equation with one spatial dimension and 
a generic linear kernel: 

dtu(x, t) = x(d x )u(x, t) — u 2 {x, t) + \j2au{x, t) r](x, t), 

with (rj(x, t)) = and (rj(xi,t 1 )r](x2,t2)) = <5(ii - £2)^(2:1 - x 2 ). 

As usual, the Langevin equation has to be taken with the Ito prescription. This means that we consider discrete 
time steps A and write 

Ui+^-mix) = x{dx)Ui{x) _ u 2 (x) + ^^- )Vi{x)j 
with the r.h.s. taken at time i and the following noise correlators: 

(r]i(x)) = and (77^ (xi)r} i2 (x 2 )) = -^S ili2 5(x! - x 2 ). 



14 The precise form of Ai (see JHJ) is not really important, we shall only ask it to be conformal invariant. 

15 We use the complex notations in which a two-dimensional vector x is equivalent to the complex number x\ + ix^- 
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We can then consider the time evolution of any functional F[ui(x)] of the event-by-event amplitude Ui(x). By series 
expansion and using the Langevin equation, we have 

F[u l+1 {x)} = F[ Ui (x)} 

+ j dx 6F S [ ^ ] A[ X (d x ) Ul (x) - uj{x) + y/2^(x)rH{x)] 

+ - / dx l dx 2 1 V n 



du^xi) 5ui{x 2 ) 
A[x(d Xl )ui(x!) - u 2 (xi) + \/2aUi{xi)rii{xi)] 
A[x(dx 2 )ui(x 2 ) - u 2 (x 2 ) + \/2<j Ui(x 2 ) r)i(x 2 )}. 

We want to obtain the evolution of the average amplitude (F[u(x, t)]). For that, we take the average in the previous 
equation and use the correlators of the noise term: 

(F[u i+1 (x)]) - (F[ Ui (x)]) 
A 

^/^( : W Ma - KW -"- (l)l ) 

If,, I S 2 F[u(x)] . n . . . . . 1 
+ 2 J dXldX2 \ Su( Xl )Sn{x 2 ) A2 "V^(xiK(x 2 ) ^S( Xl -x 2 ) 



where we have only kept the leading terms in A. 

One can finally take the continuum limit for the time variable and obtain 

d t (F[u(x,t)]) 



I 



, 6F[u(x,t)] r . , 2l 8 2 F\u(x,t)] 

dx ( -±L^Ji [x{dx )u(x,t) - u 2 (x,t)} + -^U-n^t) ) . 



In particular, one can consider the case F[u(x)] — u k (x), from which one recovers the infinite hierarchy of equations 

d t (u k (x,t)) 

= k (u k - 1 (x,t) X (d x )u(x,t)) - k (u k+1 (x,t))+k(k-l)cr (u*- 1 ^,*)). 

In this expression, we clearly see that the last term, responsible for fluctuations, directly comes from the noise term 
in the Langevin equation. The two remaining ones correspond to the Balitsky hierarchy i.e. linear evolution nd 
saturation corrections. 



APPENDIX C: SATURATION SCALE AT NONZERO MOMENTUM TRANSFER 



Let us show in more details how the reference scale for the saturation scale at nonzero momentum transfer moves 
from t = q 2 at large t to a soft scale k" 2 when t becomes small. From the condition for travelling-waves formation, 
we need to satisfy the requirement that the amplitude in the dilute tail (described by the linear equation) can be 
expressed as a superposition of waves. Hence, our starting point is the generic solution for projectile-target scattering 
in the BFKL limit: 

T lin (k,q)= J |L e a ^r(k,q)^(gT,q), 

where k is the momentum of the projectile dipole, q is the momentum transfer and Qt is a typical scale for the target. 
The crucial point at this stage of the analysis is that, in this momentum space representation, we have a factorisation 
between a target and a projectile contribution. Since we are interested in the dilute tail of the scattering amplitude, 
the scale of the projectile is bigger than Q 2 , t and Q\. In that limit, we know the asymptotic behaviour of the BFKL 
eigenfunctions 



r(k,q)«(£) 



21 



On the target side, two limits are possible depending on the ordering between t and Q\. Again, we can expand the 
BFKL eigenfunctions to get 16 

<%(<3T,q) oc 1 

Qt>9 ( Qt 



OC 



Q 2 



Reinserting everything back into the expression for the amplitude, we obtain the required exponential behaviour in 
the two limits: 



Tn„(k, q) oc / *L e *«-rW with L = * * >> 

J 2i7r ^log(/r/Qr) if q < 



This means that the reference scale for the saturation scale is fixed by the hardest between the target soft scale Q\ 
and the momentum transfer t as anticipated from the beginning. 
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